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ABSTRACT 

We present comprehensive seeing statistics for the San Pedro Martir site derived from 
the Thirty Meter Telescope site selection data. The observations were obtained be- 
tween 2004 and 2008 with a Differential Image Motion Monitor (DIMM) and a Multi 
Aperture Scintillation Sensor (MASS) combined instrument (MASS -DIMM). The pa- 
rameters that are statistically analised here are: whole atmosphere seeing measured 
by the DIMM; free atmosphere seeing -measured by the MASS-; and ground-layer see- 
ing (GL) -difference between the total and free-atmosphere seeing-. We made a careful 
data coverage study along with statistical distributions of simultaneous MASS - DIMM 
seeing measurements, in order to investigate the nightly, monthly, seasonal, annual 
and global behaviour, as well as possible hourly seeing trends. Although this cam- 
paign covers five years, the sampling is uneven, being 2006 and 2007 the best sampled 
years in terms of seasonal coverage. The overall results yield a median seeing of 0.78 
(DIMM), 0.37 (MASS) and 0.59arcsec (GL). The strongest contribution to the whole 
atmosphere seeing comes, therefore, from a strong ground layer. We find that the best 
season is summer, while the worst one is winter, in accordance with previous studies. It 
is worth noting that the best yearly results are correlated with the best sampled years. 
The hourly analysis shows that there is no statistically significant tendency of seeing 
degradation towards dawn. The seeing values are slightly larger than those reported 
before. This may be caused by climate changes. 

Key words: atmospheric effects - site testing 



1 INTRODUCTION 

In the early seventies a new observing site began operations 
at the Sierra de San Pedro Martir (SPM), Baja California, 



* Based on observations obtained at the Observatorio As- 
tronomico Nacional at San Pedro Martir, Baja California, Mexico, 
operated by the Instituto de Astronomia, Universidad Nacional 
Autonoma de Mexico, 
f E-mail: leonardo@astro.unam.mx 



Mexico. The site was selected through satellite photographs, 
and was found to be one of the three best cloud-free areas 
in the world. 

We now know that this site is one of the best astronom- 
ical locations in the world and has been considered by the 
international astronomical community as a potential place 
for large telescopes to be built in the near future being a 
candidate site of projects such as the Large Synoptic Survey 
Telescope (LSST) and the Thirty Meter Telescope (TMT), 
as well as other astronomical projects. 
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Several climatological properties have been reported, 
mainly during the first years of operation of t he Obser- 
vator i o Ast r onomico Nacional (P A N) at SPM jMendoza 
197l| . I1973I: iMendoza et al] I1972I: iTlvarez fc Maisterrena 



1 200a) 



Carrasco et al 



mo: 



19771: lAlvarez fc Loped Il982l: IWalkej |1984|) and later by 
Taoial lll99l l2003h: lEchevarna et aiT (Il998h: iHiriart et al. 
I200l| ) : iMichel et all ll2003dl ah: ICarrasco fc Sarazin 



iTapia et all (l2007ah: Alvarez et all hoOlh: iBohieas et all 
ll2008h: lOtarola et al.l 1I2OO9I): iBohigas fc NufiezI ||2010|); 



Avila et al. 



ll2006h 



lAraiza Quiiano fc Cruz-Gonzalez ( 2011 ) and Carrasco et alj 
l|2012t ). 

Results on the seeing and optical tu r bulenc e 
above the ground can be found in lAvila et al.l (Il998h: 
lEcheyarrfa et all (Il998l'): IConan et"all (|2002D :1 Avila et alJ 
l|2007h : ISanchez et all (|2007^ and 'Avila et al. ' ('2OIII), as well 



as atmosphere model l ing by [Ma sciadri & Garfia s (l200l[): 
iMasciadri et al.l (|2002l . |2003| , 120041 ) and IVogiatzis fc HiriartI 
1 20041 ) . E xtinction and opa city studies have also been 



made bv ISchuster fc Parrad (|200ll ): iHiriarj l|2003l) and 
iParrao fc Schuster! 1 20031) . Site prospec tion s tu dies within 
the SPM site have been made by ISohnI (|2007l ) and 
iBohigas et al.l (|2008l ). Comprehensive rev i ews o n th e 
site can be found in ICruz-Gonzalez et al.l l|2003l . l2004h ; 
iTapia et all l|2007al lbh and I Wehingeil l|2007h . 

In this article, we analyse the data collected by the TMT 
Project Site Survey at the SPM site and discuss and compare 
results of Differential Image Motion Monitor (DIMM) and a 
Multi Aperture Scintillation Sensor (MASS) with previous 
studies. The observations include data in the time period 
2004 O ctober to 2008 February analysed bv lSkidmore et al] 
l|2009l ). plus data taken between 2008 February to 2008 Au- 
gust whose analysis has not been presented anywhere in the 
astronomical literature, improving considerably the analysis 
of the complete TMT MASS -DIMM seeing survey at SPM 
and supplementing previous seeing results published in the 
literature. 

We also present in a detailed manner the number of 
observations and the percentage of time covered by them. 
This a novel way of treating seeing data. 

In the spirit of providing the astronomical community 
with detailed information of the seeing behaviour in SPM, 
we present nightly, monthly, seasonal, annual and global 
statistics. We too, present an analysis of the hourly see- 
ing trend which yields different results from those found by 
ISkidmore et~ai] l|2009l ). 

A summary of prior SPM seeing studies are mentioned 
in Section [2l SPM TMT site testing instrumentation and 
data coverage are presented in Section |31 followed by the 
results and detailed statistics in Section |4] 



2 THE SPM SITE: SEEING STUDIES 

The atmospheric turbulence is usually studied through a 
number of parameters, one of which is known as seeing (e). 



The relation between e, Fried parameter ro, and the turbu- 
lence integral is given by 

A 



£ = 0.1 



5.25A 



-1/5 



Cl{h)dhf', (1) 



where A is the wavelengt h and fj^ C% ( h) dh is the optical 
turbulence energy profile (|Roddier|ll98ll ). 

We note that the seeing units are given in arcseconds 
(arcsec) and since seeing is a wavelength-dependent turbu- 
lence parameter it is usually calculated for 0.5 /xm. The see- 
ing value is also corrected for the direction of observation 
and is referred to zero zenith angle. 

Early works on local SPM s e eing c onditions were done 
bv lMendozal |l973) and IWalkeil (Il97ll). and several recent 
studies have been carri ed out by lEchevarrfa et al.l (| 19981 ): 
IConan et ahl (|2002l ) and iMichel et al] ll2003al). Optical tur- 
bulence stu dies have been made b y Avila et alJ (y.998, 2004|. 
I2006l . l2007l ) and I Avila et^ l|201ll ). A co mprehensiv e revie w 
of some of these studies can be found in lEchevarn'al (|2003l ). 

The lEchevarrfa et al.l (| 19981 ) observations were obtained 
with two seeing monitors and a Micro Temperature Array 
tower (MTA); the Site Testing Telescope (STT) from Stew- 
ard Observatory was designed to observe Polaris, while the 
Carnegie Monitor (CM) could observe and track any star. 
The MTA consisted of platinum detectors, which measured 
temperature differences located at different heights, from 4 
to 28 m. These authors report a median seeing of 0.61 arcsec 
with a first quartile of 0.50 arcsec and a decrease of 0.1 arcsec 
at a height of 15 m. The STT observations covered 386 
nights, while the CM observations covered 114 nights, span- 
ning a three year period. A summary of their results is shown 
in Table [J 

IConan eFZI (l20o3 ) measured the wavefront outer scale 
and included ground-based seeing measurements with a 
DIMM monitor. They reported observations during 31 
nights. The telescope was located at two sites: at the CM 
tower and at a low altitude location. In their article they 
found a bimodal distribution during their 2000 December 
campaign with peaks centred at 0.50 and 0.75 arcsec, and an 
overall seeing median of 0.77 arcsec including all observations 
duri ng the 31-night perio d. 

IMichel et al.l (|2003al ) conducted a study with the same 
DIMM monitor during 123 nights spanning almost a three 
year period. The monitor was located at the CM tower 
at a height of 8.3 m. These authors reported a median of 
0.62 arcsec and a fir st quar tile of 0.49 arcsec. 

ISkidmore et al.l (|2009l ) analyse data collected by the 
TMT Project Site Survey at five different sites which include 
SPM. A partial set (2004 Oct. to 2008 Feb.), of combined 
DIMM and MASS data, yields seeing median values for the 
ground layer 0.58 arcsec, for the free atmosphere 0.37 arcsec 

and for the whole atmosphere 0.79 arcsec. 

Concerning optical turbulence profiling, lAvila et alj 

(199^, I2OO4I ) monitored the vertical distribution of the opti- 
cal turbulence strength using a Generalized SciDAR (GS). 
Two observation campaigns have been carried out at the 
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Table 1. SPM seeing monitoring results (whole atmosphere seeing). 



Method 


1st Quart. 

(arcsec) 


Median 
(arcsec) 


3rd Quart. 

(arcsec) 


Nights 


grprpl 


0.50 


0.61 




386 




0.48 


0.63 




99 


CM / STTi''" 


0.46 


0.58 




57 


DIMM^ 




0.50; 0.75'' 




14 


DIMM2 


0.61 


0.77 


0.99 


31 


DIMM^ 


0.48 


0.60 


0.81 


123 


MASS-DIMM* 


0.61 


0.79 


1.12 


C 


SciDARS 


0.50 


0.68 


0.97 


27 



^ lEchevarria et al.l (119981) . ^ IConan et al] | |2002|) . ^ iMichel et al] ll2003al 'l. "^ ISkidmore et al.l ||2009|) . ^ lAvila et all (|2oT3) 

° CM simultaneous with STT, Bimodal distribution, " (2004 Oct. - 2008 Feb.) 



OAN-SPM in 1997 and 2000. In 1997, the GS was installed 
at the 1.5m and 2.1m telescopes (SPM1.5 and SPM2.1) for 
8 and 3 nights (March and April), respectively, whereas in 
2000 the instrument was installed for 9 and 7 nights (May) 
at S PM1.5 and S PM2. 1, respectively. 

lAvila et al.l (|201ll ) performed a recalibration of the 
C%{h) profiles obtained in the OAN-SPM, following the 
results of a theoretical work on the nor malization proce- 
dure used in the GS data reduction method (jAvila fc Cuevaa 
20091). From a statistical analysis of the recalibrated profiles, 
Avila et al.l ([2011,) have found that the seeing in the first 
two kilometers, in the free atmosphere, and in the whole at- 
mosphere had median values of 0.38, 0.35 and 0.68 arcsec, 
respectively. 

A synthesis of the non-differential monitors (STT & 
CM), differential (DIMM) and SciDAR seeing monitoring 
campaigns is shown in Table [T] 



3 TMT SITE TESTING DATA 

3.1 The 2004-2008 TMT site testing campaign 

SPM was one of the fi ve candidate sites selected by the 
TMT site testing team (jSchock et al.ll2009l ). Three south- 
ern hemisphere sites were studied, Cerros Tolar, Armazones 
and Tolonchar in northern Chile, and two northern hemi- 
sphere sites, the 13 North (13 N) site on Mauna Kea, Hawaii 
in the United States and SPM in Mexico. During a period of 
approximately five years, from 2004 to 2008, the TMT group 
measured the atmospheric properties of each site with the 
same instrumentation and at least 2.5 annual cycles of data 
were acquir ed on each of t he ca ndidate sites. As is described 
in detail bv lSchock efal] (|2009l ) a suite of eight instruments 
were deployed in the candidate sites. The data acquisition 
methodology is also presented in their work. The acquired 
data of all the instruments is kindly available to the public 
at |http: //sitedata. tmt . org/| by the TMT organisation. 



30000 

TO 

TO 25000 
Q 

Z 20000 
c 

J 15000 

o 10000 

H 

5000 




Figure 1. Total monthly data for the 2004-2008 period. Pat- 
terns: DIMM (dense), MASS (sparse), simultaneous MASS- 
DIMM (dark). 

3.2 Seeing Instrumentation 

The SPM TMT T4 site testing station was equipped with a 
35 cm telescope mounted approximately 7 m above ground. 
At the Cassegrain focus of the telescope a combined instru- 
ment was deployed, which included a DIM M and a MASS, 
hereaf ter called MASS -DIMM described bv lKornilov et all 
(|2007l ). Combined MASS -DIMM instruments are exten- 
sively used fo r seeing measurement s and optical turbulence 
profiling (e.g. iTokovinin et al.|[20o3 . |2005| ). 

The DIMM channel measures the wavefront slope differ- 
ences over two small pupils some distance apart. This yields 
the whole atmosphere seeing from the telescop e level to the 
top of the atmosphere IjSarazin fc Roddiei<1990l ). DIMM see- 
ing measurements edimm are affected by different sources of 
bias such as image threshold, defocus, exposure time, photon 
noise, and high-frequency vibrations that must be monitored 
in real time, thus part of the measurements are rejected. In 



T 1 1 1 1 1 1 1 1 1 1 r 
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Table 2. Amount of data collected and coverage percentage for DIMM. 



Year 


Month 


2004 


2005 


2006 


2007 


2008 


Total 


January 


- 


- 


129 


0.9% 


3371 


23.1% 


6263 


43.1% 





0.0% 


9763 


16.8% 


February 


- 


- 





0.0% 


3082 


23.5% 


7767 


60.9 % 


4449 


35.5% 


15298 


30.0% 


March 


- 


- 





0.0% 


5740 


46.3 % 


5992 


47.0 % 


4794 


38.9% 


16526 


33.0% 


Winter 






















41587 


26.6% 


April 


- 


- 


2846 


25.5 % 


8185 


74.4 % 


3901 


36.4% 


9557 


86.7% 


24489 


55.7% 


iviciy 






6755 


64.9 % 


7640 


73.5 % 


8472 


81.5 % 


7657 


73.6 % 


30524 


73.4 % 


June 






5740 


59.9 % 


6609 


69.0 % 


3959 


41.4% 


6233 


65.1% 


22541 


58.8 % 


Spring 






















77554 


62.6% 


July 






628 


6.3% 


5064 


50.3 % 


5792 


57.7% 


4155 


41.2% 


15639 


38.9 % 


August 









0.0% 


6578 


60.3 % 


6174 


56.1% 


642 


6.1% 


13394 


30.6 % 


September 









0.0% 


8361 


71.2% 





0.0% 






8361 


23.8 % 


Summer 






















37394 


31.1% 


October 


6081 


45.5 % 


445 


3.2% 


8712 


66.1% 


8612 


63.7% 






23850 


44.6% 


November 


3887 


28.0% 


10134 


73.6 % 


9916 


72.1% 


8004 


57.7% 






31941 


57.8% 


December 





0.0% 


5165 


35.1% 


5933 


40.3 % 





0.0% 






11098 


18.8 % 


Autumn 






















66889 


40.4% 


Total 


9968 


36.7% 


31842 


22.4 % 


79191 


55.8 % 


64936 


45.4% 


37487 


43.4% 


223424 


40.2 % 



Table 3. Amount of data collected and coverage percentage for MASS. 



Year 


Month 


2004 


2005 


2006 


2007 


2008 


Total 


January 






241 


1.7% 


2663 


18.2 % 


5768 


39.7% 





0.0% 


8672 


14.9% 


February 









0.0% 


2448 


18.7% 


6451 


50.6 % 


2536 


20.3% 


11435 


22.4 % 


March 









0.0% 


4669 


37.7% 


4890 


38.4% 


354 


2.7% 


9913 


19.7% 


Winter 






















30020 


19.0% 


April 






2414 


21.6% 


7283 


66.3 % 


3562 


33.2% 





0.0% 


13259 


30.3 % 


May 






6164 


59.1% 


6846 


65.9 % 


7366 


70.9 % 


6106 


58.7% 


26482 


63.6% 


June 






5006 


52.3 % 


5516 


57.5% 


3266 


34.2% 


3274 


34.2% 


17062 


44.5 % 


Spring 






















56803 


46.1% 


July 






463 


4.6% 


4509 


44.8 % 


2294 


22.8% 





0.0% 


7266 


18.1% 


August 









0.0% 


5816 


53.3 % 


5837 


53.6 % 





0.0% 


11653 


26.7% 


September 









0.0% 


7641 


65.2% 





0.0% 






7641 


21.7% 


Summer 






















26560 


22.2% 


October 


5159 


38.9% 


336 


2.5% 


7223 


54.7% 





0.0% 






12718 


24.0% 


November 


3436 


25.5% 


7428 


53.7% 


6786 


49.3 % 


263 


1.9% 






17913 


32.6% 


December 





0.0% 


3908 


26.6 % 


4687 


31.8% 





0.0% 






8595 


14.6 % 


Autumn 






















39226 


23.7% 


Total 


8595 


21.5% 


25960 


18.5 % 


66087 


46.9 % 


39697 


28.8 % 


12270 


14.5% 


152609 


27.9% 
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Table 4. Amount of data collected and coverage percentage for simultaneous MASS — DIMM. 



Year 


Month 


2004 


2005 


2006 


2007 


2008 


Total 


January 


- 


- 


63 


0.4% 


2613 


17.9 % 


5451 


37.5% 





0.0% 


8127 


14.0% 


February 


- 


- 





0.0% 


2394 


18.2% 


6152 


48.3% 


2490 


20.0% 


11036 


21.6% 


March 


- 


- 





0.0% 


4521 


36.5% 


4713 


37.0 % 


352 


2.7% 


9586 


19.1% 


Winter 






















28749 


18.2% 


April 


- 


- 


2399 


21.5% 


7129 


64.9 % 


3422 


31.9% 





0.0% 


12950 


29.6 % 


iviciy 






5779 


55.4 % 


6707 


64.5 % 


7156 


68.9 % 


5946 


57.2 % 


25588 


61.5 % 


June 






4857 


50.7% 


5331 


55.6 % 


3088 


32.3% 


3117 


32.5% 


16393 


42.8 % 


Spring 






















54931 


44.6% 


July 






340 


3.4% 


4134 


41.0% 


1751 


17.6% 





0.0% 


6225 


15.5% 


August 









0.0% 


5625 


51.5% 


4511 


41.1% 





0.0% 


10136 


23.2 % 


September 









0.0% 


7360 


62.8 % 





0.0% 






7360 


20.9 % 


Summer 






















23721 


19.9% 


October 


4033 


30.5 % 


337 


2.5% 


6941 


52.6 % 





0.0% 






11311 


21.4% 


November 


3259 


24.2 % 


7274 


52.6 % 


6607 


48.0 % 


241 


1.7% 






17381 


31.6% 


December 





0.0% 


3819 


25.9 % 


4451 


30.2% 





0.0% 






8270 


14.0% 


Autumn 






















36962 


22.3% 


Total 


7292 


18.2% 


24868 


17.7% 


63813 


45.3 % 


36485 


26.4% 


11905 


14.1% 


144363 


26.4% 



particular, for the system used here, the pr ecision reported 
is better than 0.02 arcsec (|Wang et al.ll2007h . 

The MASS channel detects rapid variations of light in- 
tensity in four concentric apertures using photomultipliers 
and reconstructs a tur bulence p rofile at six different heights 
above the telescope, (jTokovinin fc Kornilovii2007i ). MASS is 
not sensitive to turbulence near the ground as this does not 
produce any scintillation and it can only measure the see- 
ing in the free atmosphere (FA), thus excluding the ground 
layer, so we have emass = £fa- From the turbulence profile 
the seeing integrated from 500 m above the telescope to the 
top of the atmosphere is computed. MAS S seeing accuracy 
is better than 0.05 arcsec (|E1s et al.ll2008l '). 



3.3 SPM TMT Data Coverage 

The SPM TMT T4 station acquired MASS -DIMM data 
from 2004 October to 2008 August. The DIMM and MASS 
nightly operations commenced one hour after sunset and 
ceased one hour before sunrise. DIMM and MASS measure- 
ments were triggered simultaneously, DIMM took integrated 
stellar light measurements during a time interval of 36 s, 
whereas MASS sampled light during 60 s periods. As the 
DIMM channel was used to acquire the target star, this re- 
sulted in s imultan eous MASS -DIMM measurements every 
70 to 90s (|E1s et al. 2009.). Sometimes computations failed 
to provide valid results probably due to cloudy conditions 



or technical problems, while sometimes the telescope system 
was shutdown due to bad weather conditions. 

In Fig.[T]we present the total amount of DIMM, MASS 
and simultaneous MASS - DIMM data for each month dur- 
ing the acquisition period. By simultaneous MASS -DIMM 
data we mean measurements within the same minute inter- 
val. 

We have also computed the statistics of the data cov- 
erage per night and per month. To compute the monthly 
temporal data coverage we have considered the duration for 
each night of the year due to the diurnal cycle variation. 
This yields a percentage value of the available data for each 
night in a given month. A 100 % value would indicate data 
obtained during the whole night for every night of the month. 

To our knowledge, this is the first time that a complete 
data coverage of seeing observations is presented in the as- 
tronomical literature. Whereas in other papers it is spoken of 
number of nights observed regardless of whether one seeing- 
point or many have been observed per night, we emphasise 
the effective length of time observed per night; making this 
a novel way of treating seeing data. 

In Tables [21 E] and [3] we give the percentage of monthly, 
seasonal and annual data coverage of DIMM, MASS and 
simultaneous MASS -DIMM during the acquisition period. 

In Figs. [2l [3] and [4] we present grey level plots showing 
only the monthly percentage of data coverage for the 2004 - 
2008 period for DIMM, MASS and simultaneous MASS- 
DIMM, extracted from Tables [1|3] and g] 
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DIMM Data Coverage (%) MASS Data Coverage (%) 




Month Month 



Figure 2. Grey level plot showing the monthly percentage of data 
coverage for the 2004-2008 period for DIMM. A black square 
means that no usable data were available. 



Figure 3. Grey level plot showing the monthly percentage of 
data coverage for the 2004-2008 period for MASS. A black square 
means that no usable data were available. 



These figures are intended to be a visual aid in order to 
see, at a glance, which months are those with the best cov- 
erage. The degree of darkening corresponds to the monthly 
percentage values presented in Tables 2, 3 and 4. 

At a first glance we immediately notice two main fea- 
tures. First, that the best coverage was obtained during 2006 
and 2007; secondly, that April, May and June are the best 
covered months. 

This is further supported by a careful analysis in Ta- 
bles [IE] and H For 2006 and 2007 the DIMM data cov- 
erage is 55.8 and 45.4% respectively, while the MASS cov- 
erage for these years is 46.9 and 28.8 % respectively. The 
simultaneous MASS -DIMM coverage is 45.3 and 26.4% re- 
spectively. The overall data coverage for all years is signifi- 
cantly less (DIMM 40.2 %; MASS 27.9 %; and simultaneous 
MASS -DIMM 26.4%). Ahhough these coverage values are 
low they reflect the enormous difficulties in obtaining see- 
ing data, even in experiments designed to ensure high data 
collec tion efficiencies, as is the case of the T MT site testing 
data (jRiddle et al.ll2006l : ISchock et al.ll201Gf l. Note that the 
best sampled year for simultaneous measurements is 2006 
(45.3 %) and the worst sampled one is 2008 (14.1 %). 

Seasonal MASS -DIMM coverage in winter is 18.2%, 
spring 44.6%, summer 19.9%, and autumn 22.3%. Seasons 
have been defined at the beginning of the month: winter 
(January, February and March); spring (April, May and 
June); summer (July, August and September) and autumn 
(October, November and December). 

In the monthly case we can see that the data coverage 
is worse in January (14.0%), December (14.0%) and July 
(15.5%), while the best sampled months are May (61.5%), 
June (42.8%) and November (31.6%). 



Simultaneous MASS-DIMM Data Coverage (%) 




Figure 4. Grey level plot showing the monthly percentage of 
data coverage for the 2004 - 2008 period for simultaneous MASS - 
DIMM. A black square means that no usable data were available. 



4 RESULTS 

We have calculated the seeing from MASS, DIMM and si- 
multaneous MASS-DIMM obtaining nightly, monthly, sea- 
sonal, annual and global statistics (first quartile, median and 
third quartile, via cumulative distributions). The behaviour 
of the measured seeing expressed as log(e) was verified to 
follow normal statistics, obtaining an expected correlation 
between seeing coverage percentage and deviations from a 
log-normal distribution. The seeing distribution is expected 
to be log-normal becau se £ is a randorn variab le that ranges 
from to infinity, and iFried fc Mever j (Il974l ') confirmed To 
as a log-normally distributed random variable. 
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Figure 5. Example of DIMM (upper curve) and MASS (lower 
curve) seeing median values during a good steady night (2007-08- 
12). 
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Figure 6. Example of DIMM (upper curve) and MASS (lower 
curve) seeing median values during a night in which only DIMM 
seeing degrades along the night (2006-03-21). 



4.1 Nightly statistics 

We have made a nightly analysis of all available data (~800 
nights of DIMM data, ~620 nights of MASS data and ~600 
nights of simultaneous data, out of ~1400 nights which cor- 
responds to the total length of the campaign), which are 
either partially covered or fully covered nights. We found a 
great variety of both DIMM and MASS behaviours; nights 
with excellent seeing throughout the night; erratic seeing 
nights; degrading seeing in which there is a sudden burst of 
bad seeing; etc. In the following paragraphs we will discuss 
a couple of examples. 

In Fig. Owe present a night (2007 August 12) with good 
steady seeing. For both DIMM and MASS the seeing vari- 
ations are well correlated. The median seeing measured by 
DIMM and MASS is 0.53 and 0.19arcsec, respectively. 

In Fig. [6] we present a night (2006 March 21) in which it 
is clear that there is a strong lack of correlation between the 
DIMM and MASS seeing. Only the DIMM seeing degraded 
during the night, while the MASS seeing remains steady with 
a median of 0.24 arcsec. 



The median standard variation of the whole atmosphere 
seeing within a night is 0.19 arcsec (with first and third quar- 
tiles 0.12 and 0.29 arcsec, respectively). 

As an example of simultaneous MASS - DIMM nightly 
statistics we show in the left panel of Figs. [7] and |8] the seeing 
median values for some of the best sampled months (2006 
September with 25 nights, and 2007 August with 24 nights). 



4.2 Monthly statistics 

In Table [S] we present detailed monthly statistics for DIMM 
and MASS data for the whole campaign. 

Apart from the worst median seeing value obtained in 
2005 January, which happens to be the month with the poor- 
est coverage (0.4%, c.f. Table |4|, DIMM seeing values vary 
from 0.6 to 1.1 arcsec and MASS seeing values vary from 0.3 
to 0.4 arcsec. An example of monthly statistics is shown in 
the right panels of Figs. [7] and [8] containing the seeing Cu- 
mulative Distribution Functions (CDF) for the same months 
previously mentioned. For DIMM we obtain a median seeing 
of 0.62 and 0.61 arcsec, respectively; while for MASS we get 
0.31 and 0.27 arcsec, respectively. 



4.3 Seasonal statistics 

We show in Fig. [9] the seeing seasonal behaviour. For each 
season we plot the DIMM and MASS median seeing values 
for the corresponding months. The DIMM seeing displays 
a pronounced seasonal variability and there seems to be a 
slight correlation between the DIMM and MASS variations. 
It is clear that the worst whole atmosphere seeing occurs in 
winter, although at the end of the season there is a tendency 
towards a better seeing. This tendency continues until the 
middle part of summer when the best seeing is achieved, 
worsening progressively towards autumn. 

In the upper part of Table [5] we present the sea- 
sonal statistics for the entire observing run. As mentioned 
above the worst whole atmosphere seeing occurs in winter 
(0.97 arcsec) and the best in summer (0.64 arcsec). 

Differences in seasonal MASS are smaller indicating that 
most of the large DIMM seeing variations are due to the 
turbulence occurring near the ground (see Sect. I4.5']l . 



4.4 Annual statistics 

We present in Figs. [10] and [11] the DIMM and MASS seeing 
measurements for the years 2006 and 2007 which are the 
best sampled years in terms of seasonal coverage. 

In both figures we can see a clear seasonal variation of 
the DIMM seeing. Whereas the MASS values tend to remain 
more stable throughout the year. 

The lower part of Table [6] presents the yearly DIMM 
and MASS statistics. It is interesting to note that the best 
whole atmosphere seeing occurs for 2006 (0.74 arcsec) which 
is the best sampled year. 
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Night of the month Seeing (arcsec) 



Figure 7. Left panel: DIMM (circles) and MASS (squares) statistics for simultaneously obtained data. Median seeing values (with lower 
and upper limits represented by first and third quartiles respectively) for each night of the month 2006 September. Right panel: Seeing 
Cumulative Distribution Function for the same month. Left curve represents MASS values with a median of 0.32 arcsec. Right curve 
represents DIMM values with a median of 0.63 arcsec. 



Table 5. Median seeing monthly values (arcsec). 



Year 



Month 2004 2005 2006 2007 2008 Total 





DIMM 


MASS 


DIMM 


MASS 


DIMM 


MASS 


DIMM 


MASS 


DIMM 


MASS 


DIMM 


MASS 


January 






1 


47 





82 





75 





36 


1 


10 





47 






1 


01 





44 


February 














1 


07 





52 





98 





41 


0.86 


0.34 





99 





41 


March 














1 


08 





33 





82 





39 


0.96 


0.50 





92 





37 


April 









98 





40 





83 





37 





92 





45 









88 





40 


May 









77 





28 





64 





30 





76 





30 


0.89 


0.43 





75 





32 


June 









84 





30 





74 





41 





73 





39 


0.65 


0.33 





75 





36 


July 









94 





59 





72 





46 





58 





30 









68 





41 


August 

















64 





35 





61 





27 









63 





31 


September 

















62 





31 

















62 





31 


October 


1.08 


0.46 





85 





33 





75 





38 

















83 





40 


November 


0.97 


0.41 





83 





42 





83 





40 





89 





52 









85 





41 


December 









78 





38 





82 





38 

















80 





38 



4.5 Ground layer statistics 



We calculated the Ground layer (GL) seeing eoL (from 7- 
500 m) which is defined as the difTerence between the DIMM 
and MASS seeing. As we note from equation [T] the contribu- 
tions of the different layers to the seeing must be summed 



as the 5/3 power, which yields 

_ ,-_5/3 _ 5/3 n3/5 /r,N 
= GL — (EoiMM ^MASS) ' 

In the data acquisition process it might occur 
that epi MM < emas s- These occurrences correspond 
to what IrokoviniTil ((2007 ) defined as "over-shoots". 
We identified that the total number of these occur- 
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Night of the month Seeing (arcsec) 



Figure 8. Left panel: DIMM (circles) and MASS (squares) statistics for simultaneously obtained data. Median seeing values (with lower 
and upper limits represented by first and third quartiles respectively) for each night of the month 2007 August. Right panel: Seeing 
Cumulative Distribution Function for the same month. Left curve represents MASS values with a median of 0.27 arcsec. Right curve 
represents DIMM values with a median of 0.61 arcsec. 



Table 6. Seasonal and annual seeing statistics (arcsec). 







DIMM 






MASS 






Ground Layer 






Season 


1st 


Median 


3rd 


1st 


Median 


3rd 


1st 


Median 


3rd 


Coverage 




quartile 




quartile 


quartile 




quartile 


quartile 




quartile 


% 


Winter 


0.67 


0.97 


1.38 


0.26 


0.40 


0.64 


0.51 


0.73 


1.11 


18.2 


Spring 


0.61 


0.78 


1.07 


0.23 


0.35 


0.54 


0.46 


0.59 


0.81 


44.6 


Summer 


0.54 


0.64 


0.79 


0.24 


0.33 


0.48 


0.38 


0.46 


0.56 


19.9 


Autumn 


0.64 


0.83 


1.16 


0.28 


0.40 


0.58 


0.46 


0.63 


0.92 


22.3 


Year 


2004 


0.72 


1.04 


1.51 


0.31 


0.44 


0.65 


0.56 


0.81 


1.27 


18.2 


2005 


0.65 


0.82 


1.13 


0.24 


0.36 


0.54 


0.49 


0.65 


0.89 


17.7 


2006 


0.58 


0.74 


1.02 


0.25 


0.37 


0.55 


0.42 


0.54 


0.76 


45.3 


2007 


0.60 


0.80 


1.17 


0.24 


0.37 


0.58 


0.46 


0.60 


0.87 


26.4 


2008 


0.61 


0.82 


1.11 


0.26 


0.38 


0.57 


0.45 


0.62 


0.86 


14.1 



rences corresponds to ^ 1.5%. We have decided to 
reject these data noting that the amount of rejected 
data is rather small and therefore overall statistics 
of the GL seeing are not affected. 

In the GL section of Table [6] we present the seasonal 
and yearly results for the GL seeing values. While in Ta- 
ble [7] we show the monthly variation of the GL seeing for 
the entire observing run. We note that the best GL seeing 



is during July, August and September months with a value 
~0.46 arcsec. 

Again, the GL behaviour is similar to the DIMM seeing: 
best in summer, worst in winter. This is due to the fact that 
the MASS values are always nearly constant showing that 
the high-altitude turbulence at SPM is small. Therefore, the 
main contribution to the whole atmosphere seeing at SPM 
comes from a strong GL. 



© 2012 RAS, MNRAS OOO.mfTSi 



10 Sanchez, Cruz-Gonzdlez et al. 



1.4 

1.2 



-o 1-0 

(0 

^ 0-8 
to 

.E 0.6 

0) 
0) 



0.2 



Winter 



January February 



March 




August September 




October November December 



Figure 9. Median seeing values (witli first and tliird quartiles as error bars) for each montli of every season covered by tlie campaign 
(2004-2008). Symbols as in Fig. [7] left panel. 



4.6 Global statistics and hourly trend 

In the left panel of Fig. [12] we show the CDF for the DIMM, 
MASS, and derived GL seeing for all simultaneous data. Ta- 
ble [2sh£ws_thB_overa11_^^ results. 

Skidmore et al.l Hooi) argue that the DIMM seeing at 



SPM rises during the second half of the night (see their 
Fig. 4) . We decided to explore the seeing behaviour and look 
for a possible degradation along the night. In Fig.[T3]we have 
plotted the number of data as a function of the hour at which 
they were taken. It is clear that before 4 UT hours and af- 
ter 11 UT hours the number of measurements taken was 
substantially smaller than within 4 to 11 UT. This justifies 



the calculation of an overall CDF in this restricted range 
presented on the right panel of Fig. 1121 The second line of 
Table |8] shows the overall results only for the same interval. 
As it is clear, the resulting statistics does not significantly 
change when using the restricted time interval. 

In Fig. [14] we plot the DIMM and MASS hourly results 
integrated over the whole campaign. The seeing seems 
to be worse at the beginning and at the end of the 
night, however the number of data (N) which sup- 
ports this assertion is rather small and it is clear 
from the figure that these data correspond to the 
longest nights of the year (i.e. winter nights where 
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Table 7. Ground layer monthly median seeing (arcsec). 




Year 



1 2 3 4 5 6 7 8 9 101112 
Month of the year 



Figure 10. DIMM (circles) and MASS (squares) simultaneously 
obtained data. Median seeing values (with first and third quartiles 
as error bars) for each month of 2006. 
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Figure 11. DIMM (circles) and MASS (squares) simultaneously 
obtained data. Median seeing values (with first and third quartiles 
as error bars) for each month of 2007. 



Month 


2004 


2005 


2006 


2007 


2008 


Total 


January 




0.97 


0.60 


0.82 




0.76 


February 






0.77 


0.76 


0.72 


0.75 


March 






0.89 


0.61 


0.67 


0.69 


April 




0.78 


0.64 


0.65 




0.66 


May 




0.65 


0.50 


0.58 


0.65 


0.58 


June 




0.68 


0.53 


0.53 


0.49 


0.57 


July 




0.64 


0.45 


0.42 




0.45 


August 






0.46 


0.48 




0.47 


September 






0.47 






0.47 


October 


0.81 


0.72 


0.56 






0.63 


November 


0.80 


0.61 


0.62 


0.61 




0.64 


December 




0.58 


0.62 






0.60 



seeing is usually worse). The points on the plot have 
error bars corresponding to 

We have made a weighted fifth grade polyno- 
mial fit to the points. It is clear that in the interval 
delimited by the dot-dot lines (corresponding to the 
shortest night of the year) the fit is essentially linear 
and horizontal, which might indicate that the seeing 
remains constant throughout the night. We know that 
in giving the data points an error bar equal to those 
points within the 4 to 11 UT interval, where N is sufficiently 
large, have a determinant influence on the results of the fit. 
However, it is correct to do it this way since we have shown 
that the statistical results in the 4 to 11 UT interval are the 
same as those for the whole interval (see Table [8]). 

In order to explore with precision the seeing be- 
haviour at the beginning and the end of the night 
we analyse the data in a different manner: integrat- 
ing with respect to the astronomical dusk and to the 
astronomical dawn, calculated for each night of the 
year. We do this since the beginning and the end of 
the night occur at a different time (UT) throughout 
the year. 

We include Figs. [151 and [161 which show the me- 
dian seeing versus hours after astronomical dusk and 
before astronomical dawn. 

Performing a linear fit to the DIMM points con- 
tained within the length of the shortest night we find 
it to be essentially flat with slopes 9.0311 x 10^* ± 
0.0016 for Fig. [15] and 0.0024 ± 0.0016 for Fig. [TBI 

The MASS points do not seem to show any 
significant tendency in either case with slopes of 
3.9400 X lO""* ± 0.0016 and -2.4118 x 10"'' ± 0.0016. 

We, therefore, conclude that the seeing at SPM 
does not statistically change either shortly after dusk 
or towards the last few hours of the nig ht, in con- 
trast with the results reported by Skidmore et ah! 
(|2009h . and in good agreement with the findings 
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Seeing (arcsec) Seeing (arcsec) 



Figure 12. MASS, DIMM and GL statistics for simultaneously obtained data. Left panel: Cumulative Distribution Functions 
(CDFs) for the whole campaign (2004-2008) including all data points. Median(MASS) = 0.37 arcsec, Median(DIMM) = 0.78 arcsec, 
Median(GL) = 0.59 arcsec. Right panel: CDFs for the whole campaign (2004-2008) including only 4-11 UT data points. Me- 
dian(MASS) = 0.37 arcsec, Median(DIMM) = 0.77 arcsec, Median(GL) = 0.58 arcsec. Left curves represents MASS values, middle curves 
represents GL values and right curves represents DIMM values. 



Table 8. Global seeing statistics (arcsec). 







DIMM 






MASS 






Ground Layer 




Data 


1st 


Median 


3rd 


1st 


Median 


3rd 


1st 


Median 


3rd 




quartile 




quartile 


quartile 




quartile 


quartile 




quartile 


All 


0.60 


0.78 


1.11 


0.25 


0.37 


0.56 


0.45 


0.59 


0.84 


4-11 UT 


0.60 


0.78 


1.10 


0.24 


0.37 


0.56 


0.45 


0.58 


0.83 



of 'Eche varn'a et al.l ('1998') and lSanchez et al.l (120031 ) 
where no clear tendency was found. 



4.7 Comparison with previous results 

Our seasonal r esults (see Sect. 14.311 ma y be compared with 
the findings of lEchevarrfa et al.l iwM) wit h the STT and 
CM monitors and by 'Mic hel et al.l l|2003bl ') with a DIMM 
monitor. The former found a summer median of 0.58 arc- 
sec while the latter found a median of 0.55 arcsec, compared 
with 0.64 arcsec obtained here. Although, we find higher val- 
ues in all seasons in this paper, summer is still shown to be 
the best season in all works. We point out that the coverage 
percentage in summer, autumn and winter, presented in this 
paper are particularly low (^ 23%). Still the median value 



fo r spring is much h igher than previous work s; 0.60 arcsec 



Echevarn'a et al.l (|l998l ) and 0.61 arcsec in I Michel et al.l 



l|2003bl ) and 0.78 (this paper). The results for autumn are 
0.68, 0.63 and 0.83 arcsec respectively, while for winter are 
0.69, 0.77 and 0.97 arcsec respectively. Thus, the results are 
consistent, if not in the absolute values, in the fact that sum- 
mer yields the best values, while winter is the worst season. 
The mean values for all three campaigns yield 0.79, 0.68, 0.61 
and 0.70 arcsec for winter, spring, summer and autumn re- 
spectively, while the overall median average is 0.70 arcsec. 
This value is higher than the median seeing in the first 
two campaigns and lower than the results in this work (see 
Sect. 14. 6p . The differences show the skewness produced by 
uneven coverage during the year in all three campaigns. 

Ou r global seeing results are in good agreement with 
those of lSkidmore et all l|2009l ). The slight differences found 
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Figure 13. Total number of hourly simultaneous DIMM and 
MASS measurements (N). The vertical lines delimit the interval 
between 4 and 11 UT hours. 



could be due to a larger temporal coverage in our data set, 
from 2004 Oct. to 2008 Aug. 

Comparing the results in this work (Table [8| with those 
of previous studies (Table [T] and Sect. 14. 3p reveals that the 
whole atmosphere seeing {edimm) is larger than most of 
the values reported before. This could be due to a slight 
degradation in seeing in recent years. 

Based on the results for Paranal and La Silla 
ijSarazin et all l2008l fl we suggest that the appar- 
ent seeing degradation found at SPM might be 
correlated with climate change. We have consulted 
the Climatic Data Base of the Northwest Mexico 
(http://peac-bc.cicese.mx/datosclim) which con- 
tains weather data for a number of stations near the SPM 
observatory site for the years 1981 to 2008. For this time 
period, we found that the temperature tends to increase 
0.05 Celsius per year while the water precipitation tends to 
decrease 5 mm per year. These variations clearly point to a 
climate change which might produce local effects. Whether 
the seeing degradation is due to these effects is yet to be 
established. However, a full study of the dependence of the 
seeing value with the SPM climate change is beyond the 
scope of this paper. 

From a statistical point of view we found that the seeing 
at SPM does not suffer any degrad ation towards dawn, i n 
contrast with the results reported bv lSkidmore et al] (|2009h . 
However, the consistent increase of seeing values towards the 




UT (Hours) 



Figure 14. Hourly seeing computed from total simultaneous 
DIMM (upper curve) and MASS (lower curve) data. The error 
bars correspond to The fitted curve to the DIMM data 
corresponds to a weighted fifth degree polynomial. Ver- 
tical dot-dot lines indicate the shortest night, whereas 
vertical dash-dot-dot lines delimit the longest night of 
the year. 



end of the night might suggest, as ISkidmore et al.l (|2009l ) 
states, that seeing degrades towards dawn. To establish this 
it would be necessary to collect a large number of observa- 
tions towards dawn so that this result may be established 
wit h a sufficient degree of statistical significance. 

ICarrasco et all dioil l made a study of SPM solar ra- 
diation (cloud coverage) using the TMT data for the same 
period presented in this paper. They state that SPM skies 
are clear in spring, relatively cloudy in summer and fairly 
clear in winter. It is interesting to notice that there appears 
to be a correlation between the cloud cover during day time 
and the nightly value of the seeing, in the sense that cloudy 
days have better seeing at night and viceversa. This might 
indicate that the processes that occur when the cloud cover 
is dissipated contribute to smaller amounts of atmospheric 
turbulence. This is an interesting point which merits further 
study that is beyond the scope of this work. 



5 SUMMARY AND CONCLUSIONS 

In this work we have used the data obtained by the TMT Site 
Testing Project for the SPM site in the north-west of Mexi- 
co. The data consist in DIMM and MASS measurements 
taken over a period of nearly five years (2004 October - 
2008 August, ^1400 nights). A variety of previous studies 



have measured the seeing at SPM. The obtained median 
1 Sarazin M., 2010, |http://www.eso.org/gen-fac/pubs/astcUm/paranal/se aiiiitg^3ragtoar |Etotti^ .50 to 0.79arcsec; although these mea- 
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Figure 15. Hourly seeing computed from total simultane- 
ous DIMM (upper curve) and MASS (lower curve) data, 
for hours after astronomical dusk. The error bars corre- 
spond to Vertical dot-dot lines delimit the shortest 
night of the year. 



Figure 16. Hourly seeing computed from total simultane- 
ous DIMM (upper curve) and MASS (lower curve) data, 
for hours before astronomical dawn. The error bars cor- 
respond to Vertical dot-dot lines delimit the shortest 
night of the year. 



surements have been taken in different epochs and with dis- 
parate instruments. The coverage of each observing run is 
also limited. 

For the analysis presented in this paper, we have made 
a detailed coverage study obtaining percentages of total ob- 
served time per night with respect to total night length. 
This is the first time that this approach is presented in 
the seeing astronomical literature. The total coverage for 
DIMM, MASS and simultaneous data is 40.2%, 27.9%, and 
26.4% respectively. This means that, of the 1400 available 
nights, the DIMM data were obtained in an equivalent of 563 
"fully observed" nights (nights with a 100 % time coverage), 
whereas for the MASS this corresponds to 391 nights, while 
for the simultaneous data it corresponds to 370 nights. 

We calculated nightly, monthly, seasonal, annual and 
global statistics of the seeing from MASS, DIMM and 
simultaneous MASS -DIMM observations. The simultane- 
ous MASS -DIMM results indicate that the best see- 
ing is obtained for 2006 September: 0.62 (DIMM), 0.31 
(MASS)arcsec and 2007 August: 0.61 (DIMM), 0.27 
(MASS) arcsec. These months correspond to our definition 
of the summer season in which the calculated seeing is 0.64, 
0.33 and 0.46 arcsec, median values for DIMM, MASS and 
ground-layer respectively. Making the summer the best see- 
ing season at SPM. It is worth noting that the seeing ob- 
tained for 2006 - 0.74 (DIMM), 0.37 (MASS) and 0.54 arcsec 
(GL) - is the best yearly seeing of the whole campaign which 
also corresponds to the best sampled year. For the GL seeing 
the best month is July with a value of 0.45 arcsec. The overall 
results yield a median seeing of 0.78 (DIMM), 0.37 (MASS) 



and 0.59 arcsec (GL). Comparing with previous works we 
found that our whole atmosphere seeing values are slightly 
larger than most of the values reported before, perhaps this 
is due to a degradation of seeing with time, which might 
be caused by climate effects. The hourly analysis clearly 
showed that there is no statistically significant ten- 
dency of seeing degradation towards dawn, in con- 
trast with that reported by Skidmore (2009). 

These results show that SPM is a competitive astronom- 
ical site, in which it would be worth performing a continuous 
astroclimate evaluation. 
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